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We present a study of neutrino-nucleus interactions at the T2K experiment based on the GiBUU 
transport model. The aim of T2K is to measure v e appearance and #13, but it will also be able 
to do a precise measurement of disappearance. The former requires a good understanding of 
7T° production while the latter is closely connected with a good understanding of quasielastic scat- 
tering. For both processes we investigate the influence of nuclear effects and particular final-state 
interactions on the expected event rates taking into account the T2K detector setup. 
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I. INTRODUCTION 

After neutrino oscillations were first observed in at- 
mospheric and solar neutrino rates P, 0, an extensive 
experimental program has started aiming at the precise 
determination of v masses and mixing angles. Of the 
three mixing angles, #13 is still unknown, only upper lim- 
its have been placed so far. Still unknown is also the sign 
of one of the two mass squared differences, namely Am| 2 , 
and thus the mass hierarchy (cf., e.g., the review article 
in Ref. 0). 

The T2K experiment [|J, which has just started its 
operation, is a long-baseline high-precision neutrino- 
oscillation experiment. Its major aim is to measure v e ap- 
pearance and thus #13 but also to provide a precise mea- 
surement of disappearance to improve $23 and Ato^. 
Its target material - common to all modern detectors - 
consists of heavy nuclei. This causes a major difficulty: 
Particles produced in neutrino interactions can reinter- 
act before leaving the nucleus and can be absorbed or 
change their kinematics or even their charge before be- 
ing detected. Nuclear reinteractions limit our ability to 
identify the reaction channel by changing the topology of 
the measured hadronic final state. Consequently, the de- 
tected rates on nuclei are changed significantly compared 
to the ones on free nucleons. 

Appearance experiments like T2K search for a specific 
neutrino flavor in a neutrino beam of different flavor. The 
flavor of the neutrino can only be determined from the 
charged lepton it produces in the interaction. ir° produc- 
tion events in neutral current (NC) reactions are a source 
of background in v e appearance searches in a beam 
because they might be misidentified as charge current 
(v ei e~) interactions. A precise and well-tested model for 
NOV production on nuclei is thus necessary 

The oscillation probability depends directly on the 
neutrino energy: disappearance experiments search 
for a distortion in the neutrino flux in the detector posi- 
tioned far away from the source. From the flux difference, 
one gains information about the oscillation probability 
and with that about 623 and Am^- However, the neu- 
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trino energy cannot be measured directly but has to be 
reconstructed from the final-state particles that are de- 
tected. Usually, LBL experiments use charged current 
quasielastic (CCQE) scattering events for the energy re- 
construction. Their identification is also influenced by 
nuclear effects which have to be understood and consid- 
ered in the experimental analysis. Thus, to extract the 
oscillation parameters from the measured particle yields, 
the experimental analyses have to rely on models for the 
neutrino-nucleus interaction (see, e.g., the proceedings of 
the NUINT conference 0]). 

In the T2K experiment, JPARC's high intensity 
beam is directed toward the Super-Kamiokande water 
Cherenkov detector which is located about 300 km away. 
In addition, T2K has a high-resolution near detector 
complex to determine the energy spectrum, the flavor 
content of the beam and also neutrino cross sections in 
particular by measuring NC71- and CCQE on oxygen. 
The energy of the neutrino beam extends up to a few 
GeV with a maximum at about 600 McV. 

In T2K's case, the near and the far detectors are 
conceptually very different (water Cherenkov vs. track- 
ing detector) which implies very different event selection 
schemes. Thus, a quantitative understanding of neutrino- 
nucleus cross sections is mandatory. 

While the Giessen Boltzmann-Uehling-Uhlenbeck 
(GiBUU) transport model has already been widely ap- 
plied to general questions of neutrino interactions with 
nuclei [6j-l8( and has also looked in some more detail at 
present experiments such as MiniBooNE and K2K 0, 
in this paper we primarily investigate the implications of 
the different detector techniques on measured NC7r° and 
CCQE spectra at T2K taking into account in-medium 
modifications and final-state interactions. 



II. NEUTRINO-NUCLEUS SCATTERING IN 
THE GIBUU MODEL 

In-medium modifications and final-state interactions 
inside the target nucleus are known to be the major 
source of systematic errors in LBL neutrino experiments 
Q. Thus, the theoretical understanding of nuclear ef- 
fects is essential for the interpretation of the data and 
only by applying state-of-the-art models it is possible to 
minimize the systematic uncertainties in neutrino fluxes, 
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FIG. 1: (Color online) Neutrino energy flux for the T2K ex- 
periment. The solid line taken from Fig. 6 (b) of Ref. [lrj ] 
shows the unoscillated flux while the dashed line shows the 
oscillated flux, i.e., the flux multiplied by (1 — P OS c) (see text 
for details). 

backgrounds and detector responses. 

For this aim, we apply the GiBUU transport model 
which is a unified framework for the description of a wide 
range of nuclear reactions from heavy-ion collisions to 
pion and electron/photon scattering and has been widely 
tested and validated [HHH (see Ref. [H for more ap- 
plications). Within this approach, we treat neutrino- 
nucleus scattering as a two-step process. In the initial- 
state step, the neutrinos interact with nucleons embed- 
ded in the nuclear medium. In the final-state step, the 
produced particles are propagated through the nucleus 
undergoing complex final-state interactions (FSI). Here, 
we present only a short review of our model and refer to 
Ref. (EH! for details. 

The energy distribution of the T2K beam is shown in 
Fig. Q] by the solid line. The dashed line shows the effect 
of neutrino oscillations on the flux: The probability that 
v„ remains as ia, is 

P(v„ -> = 1 - P osc = 1 - sin 2 20 2 3 sin 2 f^f^) , 

(1) 

using the atmospheric mixing parameters 623 = 45°, 
Am 2 3 = 2.5 x 10~ 3 cV 2 , and a distance of L = 295 
km. 

In this energy region, the elementary vN reaction is 
dominated by two processes: quasielastic scattering and 
the excitation of the A resonance (P33(1232)). The vec- 
tor form factors are well constrained by electron scatter- 
ing data (see Ref. Q for details) while we have a signif- 
icant uncertainty in the axial form factors: Goldberger- 
Trciman relations have been derived for the axial cou- 
plings, but they do not give information about the Q 2 
dependence. The Q 2 dependence of the A axial form 
factor is fitted to either ANL or BNL bubble chamber 
neutrino-scattering da/dQ 2 data for the v^p —t ^n + p 
reaction. Fig. [5] shows the integrated cross section to- 
gether with the data. Note that the solid curve fits the 
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FIG. 2: (Color online) Total v^p — ¥ (i~TT + p cross section as a 
function of the neutrino energy compared to the pion produc- 
tion data of of ANL (Refs. Q3] (•), [U (■)) and BNL ([H[ 
(x)). The solid line has been obtained with a form factor 
fitted to the ANL data, the dashed line is fitted to the BNL 
data. 



ANL data (modified dipolc form factor) , while the dashed 
curve fits the BNL data (dipole form factor). Thus, the 
latter would obviously lead to higher pion production 
cross section also on the nucleus. We refer the reader 
to Chapter 5 of Ref. [H for an extended discussion. Re- 
cently, a combined form factor fit to both ANL and BNL 
data has also been performed in Ref. [13, Gil • 

In the nuclear medium, the neutrino-nucleon cross sec- 
tions are modified. Bound nucleons are treated within a 
local Thomas- Fermi approximation. They arc exposed to 
a mean-field potential that depends on density and mo- 
mentum. We also take Pauli blocking and the medium- 
modified spectral functions of the outgoing hadrons into 
account. Our model for neutrino- (bound)nucleon scat- 
tering is described in detail in Ref. Q. 

After the initial neutrino-nucleon interaction, the pro- 
duced particles propagate through and out of the nucleus. 
During propagation they undergo complex FSI simulated 
with the coupled-channel GiBUU transport model (for 
details, see Ref. [l5| and references given there). It mod- 
els the space-time evolution of a many-particle system 
in a mean-field potential including clastic and inelastic 
collisions between the particles, and also particle decays 
into other hadrons. 

To summarize, the GiBUU model provides a theory- 
based, consistent and well-tested treatment of in-mcdium 
modifications and final-state interactions. It successfully 
describes electron-and photon-induced reactions over a 
wide energy range. Because these reactions are quite 
similar to the vA — > X reaction in that the incoming 
particle interacts with all target nucleons and the vector 
couplings are the same, we consider this as an impor- 
tant benchmark for the description of neutrino-induced 
reactions. 
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FIG. 3: (Color online) Panel (a) shows the 'V e appearance 
signal" event in a Cherenkov detector while panels (b) and 
(c) are possible backgrounds. Taken from Ref. [%|. 
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The main task in a v e appearance experiment as T2K 
is to detect electron neutrinos in a (almost) pure Vp beam 
in the Super-Kamiokande detector. The signal event is 
the v e CCQE interaction as indicated in Fig. [3] (a). How- 
ever, not all of the detected v e events have their origin in 
neutrino oscillations but are already present in the initial 
neutrino beam, e.g., through muon and kaon decays. A 
good understanding of the neutrino beam composition is 
therefore essential for v e appearance searches. 

Even more significant at low reconstructed E v arc 
misidentificd events, mainly because a Cherenkov detec- 
tor like Super-Kamiokande cannot distinguish between a 
photon and an electron. Thus, induced NC ir° pro- 
duction, 1 where the 7r° decays into two 7s, is the major 
source of background when one of the photons is not 
seen or both Cherenkov rings overlap. This process is 
illustrated in Fig. [3] (b). Additional background comes 
from the excitation of a A resonance via a NC interac- 
tion followed by its radiative decay, A — > 7 AT, which also 
leads to a final state with a photon (see Fig. [3] (c)). 2 

A proper understanding of neutrino-induced pion pro- 
duction, especially NC7r°, is therefore essential in v e ap- 
pearance experiments. T2K requires the NC7r° cross sec- 
tion to be known to 10% accuracy for the resulting error 
on the oscillation parameters to be comparable to that 
from statistical uncertainties (24[. This has triggered a 
lot of experimental activity toward direct NC71- cross sec- 
tion measurements in the last years to be used as direct 
input for the oscillation analysis [2oT - l28l |. 

In the following we study the influence of nuclear ef- 
fects and final-state interactions on NC7r° production (see 
also Ref. [29|). Fig. g] shows the results for NC single- 
7T° production off 16 averaged over the incoming, un- 
oscillated T2K energy distribution Comparing the 
dashed with the solid line (results without FSI and spec- 
tral function vs. full calculation), one finds a significant 
difference. The shape is caused by the energy depen- 
dence of the pion absorption and rescattering cross sec- 
tions. Pions are mainly absorbed via the A resonance, 



1 We recall that in NC events only the outgoing hadrons and/or 
their dec ay p roducts are seen, but not the outgoing neutrino. 

2 See Ref. 1231 for a brief discussion of the radiative A decay. 



FIG. 4: (Color online) NC induced single-7r production on 
16 O as a function of the pion momentum averaged over the 
T2K flux. The dashed and the solid lines denote the calcula- 
tion with a modified dipole form factor (fitted to ANL data) 
for the A resonance. The dash-dotted and dotted lines are 
obtained with a dipole form for the form factor (fitted to the 
BNL data). 

i.e, through tt N -> A followed by AN -> NN. This ex- 
plains the reduction in the region around p w = 0.2 — 0.5 
GeV. Pion elastic scattering irN — > ttN reshuffles the pi- 
ons to lower momenta and leads also to charge exchange 
scattering into the charged pion channels. Furthermore, 
we show in this plot the result for our second form fac- 
tor set, namely the dipole form factor (dash-dotted and 
dotted lines). The difference between the two sets re- 
flects the uncertainty caused by the discrepancy in the 
old ANL and BNL data. 

The total flux averaged NCl7r° cross section is (0.34 — 
0.48) • 10 _38 cm 2 where the range again gives the uncer- 
tainty from the axial form factor. 

The vast majority of the pions come from initial A 
excitation, pion production through higher resonances 
plays almost no role. ir° can be produced in our approach 
not only through resonances but also through initial NC 
elastic scattering, i.e., final-state nuclcons can rescatter 
in the nucleus and produce pions; this accounts for w 
2%. Missing in our approach is deep-inelastic scattering 
and any non-resonant single- 7r/double-7r etc. background. 
However, we do not expect a significant contribution of 
those at the rather low T2K energy. Coherent pion pro- 
duction is also possible, but its contribution is only of the 
order of (0.02 - 0.04) • 10~ 38 cm 2 [!}. Thus, our model 
accounts for the relevant contributions. 

A crucial problem in any such experiment with 
Cherenkov counters is that the decay photons of neutral 
pions can be misidentified as electrons. We therefore, 
investigate in Fig. [5] the probability that a misidentificd 
7T° is counted as a v e appearance event. The probabil- 
ity that a 7T° cannot be distinguished from e^- is given by 
the dashed line (taken from Fig. 2 of Ref. [3l|). ^ mainly 
depends on the opening angle of the two decay photons. 
The solid line shows the weighted cross section averaged 
over the T2K flux and calculated using the elementary 
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FIG. 5: (Color online) NC induced single- 7r° production cross 
sections on 16 O averaged over the T2K flux and multiplied 
by the misidentification probability (dashed line, taken from 
Ref. [3l|]) as a function of the pion momentum. 



ANL data as input. The total cross section for misiden- 
tified events is now 0.09 • 10~ 38 cm 2 and thus 26% of the 
true pion events. 



IV. CCQE 

The other major aim of the T2K experiment is the 
precise measurement of disappearance and the cor- 
responding mixing parameters. This task requires an 
equally precise measurement of the neutrino flux both 
at the near and the far detector. From the flux differ- 
ence, one can extract the oscillation parameters. The 
neutrino energy is reconstructed using the CCQE reac- 
tion, thus, the related challenge is to identify true CCQE 
events in the detector, i.e., muons originating from an ini- 
tial QE process v^n fJ,~p. The difficulty comes from 
the fact that the true CCQE events are masked by FSI 
in a detector built from nuclei. In T2K the QE event 
selection is different in the near and far detector: in gen- 
eral, the near detector counts all events as CCQE-like 
that have a single proton track but no pions. In the wa- 
ter Cherenkov detector Supcr-Kamiokandc, serving as far 
detector, CCQE-like means, that no pion is seen in the 
event. In both detectors the FSI lead to misidentified 
events, e.g., an initial A whose decay pion is absorbed or 
undergoes "pion-less decay" looks "CCQE-like" . 

In the following, we investigate the quality of these 
event selection methods considering also pion and pro- 
ton kinetic energy detection thresholds. Fig. [5] (a) shows 
the event selection for the near detector (tracking de- 
tector: CCQE-like = proton but no pion) while the far 
detector is addressed in panel (b) (Cherenkov detector: 
CCQE-like = no pion). The flux averaged cross sections 
are plotted versus the proton kinetic energy detection 
threshold. The "true CCQE" events are denoted with 
the solid lines, the CCQE-like events by the dashed ones. 
In addition, those are separated in the two main con- 
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FIG. 6: (Color online) Total flux averaged CCQE cross section 
on le O (solid lines) versus the proton kinetic energy detection 
threshold. Shown are besides the true CCQE cross section 
two different methods on how to identify CCQE-like events 
in experiments (dashed lines). Panel (a) shows the method 
applied in the near detector (tracking detector) , panel (b) the 
one applied in the far detector (Cherenkov detector). Panel 
(c) shows the result for the far detector but averaged over 
the oscillated flux shown in Fig. Q] The contributions to the 
CCQE-like events are also classified (CCQE-like from initial 
QE (dash-dotted lines), from initial A (dotted lines)). 



tributions, namely CCQE-like from initial QE scattering 
(dash-dotted lines) and CCQE-like from initial A excita- 
tion (dotted lines). 

In the near detector considerably fewer CCQE-like 
than true CCQE events are detected (panel (a), the dif- 
ference between the dashed and the solid line) . The final- 
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FIG. 7: (Color online) Total flux averaged CCQE cross section 
on 16 O (solid lines) versus the charged pion kinetic energy 
detection threshold. The lines are as in Fig. [6] 



state interactions of the initial proton lead to secondary 
protons, or, via charge exchange to neutrons which are 
then not detected as CCQE-like any more (trigger on sin- 
gle proton track). However, the amount of fake events in 
the CCQE-like sample is small (the dashed and dash- 
dotted lines are close). Thus, the observation of the pro- 
ton helps to restrict the background, but also leads to an 
underestimate of the true CCQE cross section by at least 
25%. This difference has to be reconstructed with event 
generators that have to be very realistic in describing the 
in-medium nuclcon-nucleon interactions. The results are 
obviously very sensitive to the proton kinetic energy de- 
tection threshold because the event selection procedure 



explicitly triggers on a single proton. At a realistic de- 
tection threshold of 0.2 GeV proton kinetic energy only 
about 30% of the total QE cross section is measured and 
70% has to be reconstructed. 

The far detector is able to detect almost all true CCQE 
(panel (b), solid vs. dash-dotted lines agree approxi- 
mately) but sees also a considerable amount of "fake 
CCQE" (or "non-CCQE") events (panel (b), the dashed 
line is roughly 15% higher than the solid line). These fake 
events have to be removed from the measured event rates 
by means of event generators, if one is interested only in 
the true QE events. It is obvious that this removal is 
the better the more realistic the generator is in handling 
the in-medium it — N — A dynamics. For this detector 
the results do not depend on the proton kinetic energy 
threshold because the event selection is independent of 
the proton. 

Fig. [6] (c) shows a similar scenario as (b) but now aver- 
aged over the oscillated T2K flux. We find almost no dif- 
ference of the total CCQE-like cross section compared to 
panel (b), however, the contributions are slightly shifted 
(more A, less true QE). This is caused by the two-bump 
structure of the oscillated flux shown in Fig. [TJ The lower 
energy bump leads mainly to true CCQE events but — 
because of the lower energy — with reduced cross sec- 
tion (solid line). However, the pion cross section and 
with that the CCQE-fake cross section is increased ow- 
ing to the higher energy bump (dotted line). However, 
overall, the effect of neutrino oscillations on the CCQE 
cross section is only minor. 

Fig. [7] shows the same cross sections as the previous 
plot now versus the pion kinetic energy detection thres- 
hold. Both detector types require 'no pions' for an event 
to be CCQE-like, thus, the CCQE-like cross section in- 
creases with increasing threshold because then more and 
more pions are not seen in the detector. More exactly, 
the A contribution to the CCQE-like cross section (dot- 
ted lines), i.e., the fake events, while the QE contribution 
(dash-dotted lines) stays constant. Also here the average 
over the oscillated flux (panel (c)) is very similar to the 
one over the unoscillated flux (panel (b)), and again the 
single contributions which add up to the dashed line are 
slightly different. 

We note that the uncertainty from the A axial form 
factor is only of minor importance here because the re- 
sults are dominated by QE scattering. 

Finally, we discuss the implications for the neutrino 
energy reconstruction which is of major importance for 
Vfj, disappearance measurements. The neutrino energy is 
reconstructed from the CCQE-like events. The preced- 
ing discussion and also our detailed study in Ref. fioj 
have shown that that the QE identification is worse in 
Cherenkov detectors than in tracking detectors: The 
CCQE-like sample in the T2K far detector contains a 
significant amount of non-QE induced events while the 
CCQE-like sample in the T2K near detector is rather 
clean but significantly reduced. The neutrino energy 
is reconstructed from the measured muon properties of 
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the full CCQE-like sample assuming quasifrcc kinemat- 
ics. The non-QE induced fraction of the CCQE-like sam- 
ple causes a shift of the reconstructed energy to values 
lower than the true energy [lj| ■ This effect is relevant in 
the far detector. Because of the conceptional difference 
between the two detectors, one reconstructs different en- 
ergy spectra in the near and in the far detector with a 
larger systematic error in the latter. In Ref. (Io| we have 
shown that at a neutrino energy of 0.7 GeV, correspond- 
ing roughly to the flux- maximum in the T2K experiment, 
the standard deviation for the reconstructed energy from 
the true energy is 18% for the tracking (near) detector, 
whereas it is 23% for the Cherenkov (far) detector. Be- 
cause for an extraction of the oscillation parameters the 
flux has to be compared at the same energy, the over- 
all combined standard deviation of the energy is thus 
a/18 2 + 23 2 % = 29%. This inaccuracy enters directly 
into the determination of the neutrino oscillation param- 
eters which underlines the need for a quantitative under- 
standing of the described effects. 

V. CONCLUSIONS 

In this work we have applied the GiBUU hadronic 
transport model to the T2K experiment. For their 
physics goals, T2K needs a good description of both the 
NC7T and the CCQE cross section. The accuracy of the 
former depends strongly on the choice for the A axial 



form factor. The NC7r° spectra are strongly modified by 
FSI so that a good understanding of these is also manda- 
tory in the experimental analysis. Furthermore, our anal- 
ysis shows that about 26% of all pion events contribute 
to misidentified electron appearance events. 

We have furthermore shown that the measured CCQE 
cross section depends strongly on the detector setup and 
that it is closely entangled with the CCl7r + cross sec- 
tion. This has consequences for the neutrino energy re- 
construction in T2K where near and far detectors are 
conceptually different: The QE identification is worse 
in Cherenkov detectors than in tracking detectors and, 
therefore, the reconstructed energy has a wider error in 
the far than in the near detector. We find that an error of 
about 29% is to be expected for the reconstructed energy 
in the oscillation parameter extraction. 

To conclude, in-medium effects and, in particular, 
final-state interactions influence the measured rates sig- 
nificantly. Their understanding — within a consistent 
and well-tested model — is crucial for the data analysis 
and interpretation of the T2K experiment. 
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